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Abstract 

It  is  shown  how  permittivity  can  be  extracted  via  time  domain  reflection  data 
from  a  perfect  electric  conductor  (PEC)  backed  planar  sample  of  a  low-loss,  non- 
dispersive  dielectric  using  two  rectangular  Ku-band  waveguide  aperture  probes  with 
attached  PEC  flange  plates  of  the  same  geometry  but  different  dimensions.  Of  critical 
importance  is  being  able  to  identify  the  reflection  from  the  edge  of  the  flange  plate  in  the 
parallel  plate  region  created  by  the  plate  and  the  PEC  backing  of  the  sample.  A  signal 
processing  method  that  takes  advantage  of  physical  insight  into  the  geometry  and 
superposition  is  developed  for  identifying  this  edge  reflection. 

Measurements  are  taken  using  square  and  circular  plate  geometries.  Measured 
data  is  processed  using  both  a  Kaiser  window  and  the  signal  subtraction  method 
developed  in  this  research  to  verify  the  analysis.  Final  results  are  presented  and  future 
work  is  discussed. 
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A  Simple  Non-Destructive  Method  for  Characterizing  Non-Dispersive,  Low-Loss 

Dielectrics 

1.  Introduction 

Material  characteristics  define  how  forces  interact  with  matter.  Within  the  realm 
of  electromagnetics  (EM),  electric  and  magnetic  forces  interact  with  media  resulting  in 
conduction,  polarization  and  magnetization  currents.  These  current  densities  are  directly 
related  to  conductivity,  permittivity  and  permeability  respectively. 

The  exploitation  of  these  parameters  has  been  instrumental  in  modern  technology. 
Being  able  to  control  how  electromagnetic  fields  interact  with  a  specific  material  has 
proven  to  be  useful  in  applications  ranging  from  communications  to  cooking.  In  order  to 
make  material  parameters  useful  in  applications,  it  often  becomes  of  critical  importance 
to  ensure  that  the  constitutive  parameters  are  within  specified  tolerances.  In  order  to 
assess  the  permittivity  and  permeability  (and  corresponding  tolerances)  of  a  material,  a 
variety  of  parameter  extraction  techniques  have  been  developed  over  the  years. 

However,  a  large  portion  of  these  methods  are  either  destructive  to  the  material  under 
test,  computationally  demanding  or  both.  This  leads  to  increased  cost  since  there  is  time 
lost  in  sample  preparation,  testing  and  computation  of  the  constitutive  parameters. 
However,  methods  that  can  provide  non-destructive,  computationally  simple  ways  to 
determine  the  constructive  parameters,  even  if  limited  to  just  determining  the  permittivity 
of  a  material,  would  be  of  enormous  benefit  as  it  would  cut  costs  without  sacrificing 
quality. 
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Figure  1-1.  Basic  setup  of  probe 

1.1  Problem  Statement 

As  previously  stated,  a  majority  of  methods  for  characterizing  the  parameters  of  a 
material  are  destructive,  computationally  demanding  or  both.  In  most  cases,  in  order  to 
save  time  using  computational  resources,  a  closed  form  method  for  determining 
permittivity  and  permeability  is  used.  One  common  method  is  the  Nicholson  Ross  Weir 
(NRW)  algorithm  [2],  [6],  [16].  The  NRW  method  uses  an  exact  formulation  to  find  the 
complex  permittivity  and  permeability  of  a  material;  however,  it  is  a  destructive  method. 
This  means  that  the  material  under  test  has  to  be  machined  to  fit  in  the  waveguide  sample 
holder.  This  comes  with  two  drawbacks,  namely,  a  limited  sample  interrogation  area  and 
the  added  time  associated  with  sample  preparation. 

The  goal  of  this  thesis  is  to  develop  a  new,  computationally  efficient,  non¬ 
destructive  method  for  characterizing  a  material  using  a  simple  time  domain  technique. 

In  this  new  technique,  reflections  from  the  aperture  and  the  edges  of  a  flanged  rectangular 
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waveguide  probe  are  utilized  to  extract  the  complex  permittivity  of  a  conductor  backed 
dielectric  material  as  depicted  in  Figure  1-1. 

1.2  Limitations 

The  time-domain  analysis  technique  introduced  in  this  research  does  have  some 
limitations.  First,  since  the  edge  reflection  from  the  plates  must  be  observed,  the  material 
under  test  must  be  low  loss;  otherwise,  the  reflection  would  never  be  observed. 
Additionally,  the  material  must  have  very  low  dispersion;  otherwise,  the  temporal 
response  may  lack  sufficient  resolution  required  by  the  technique.  Furthermore,  since  the 
network  analyzer  (NWA)  does  not  collect  time-domain  data  but  uses  a  Fourier  Transform 
(FT)  to  convert  frequency-based  data  into  the  time-domain,  there  is  a  limit  in  resolution 
which  is  directly  related  to  the  bandwidth  of  the  probe  being  used  to  make  the 
measurements. 

Other  tradeoffs  are  made  by  using  only  reflection  measurements.  Since,  as  a  rule 
of  thumb,  one  can  only  solve  for  as  many  unknowns  as  measurements  taken,  it  becomes 
apparent  that  either  the  permittivity  or  permeability  can  be  computed,  but  not  both. 

1.3  Scope 

Many  methods  are  available  for  measuring  the  constitutive  parameters  of  a 
material.  However,  the  focus  of  this  thesis  will  be  on  using  waveguide  probes.  As 
mentioned  before,  the  bandwidth  limits  the  resolution  in  the  time  domain.  Therefore,  Ku 
Band  waveguides  have  been  selected  for  use.  This  decision  was  based  on  two  factors. 
First,  the  NWA  used  in  this  research  only  has  a  range  from  10.0MHz  to  20.0GHz,  and 
secondly,  within  that  range,  the  Ku  Band  provides  the  largest  bandwidth  for  a  waveguide, 
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covering  from  10.0GHz  to  18.5GHz.  A  full  development  of  how  this  bandwidth  was 
calculated  is  located  in  appendix  C.  Additionally,  since  both  constitutive  parameters 
cannot  be  simultaneously  computed,  the  scope  will  be  limited  to  non-magnetic  materials 
(i.e.,  dielectrics).  While  the  scope  has  been  limited  to  these  conditions,  the  overall 
analysis  will  be  accomplished  in  a  way  that  could  be  applied  to  other  probe  devices  with 
different  bandwidths. 

1.4  Thesis  Organization 

Chapter  2  provides  a  theoretical  background  starting  with  Maxwell’s  equations 
and  constitutive  parameters,  then  discusses  wave  behavior  in  parallel  plate  systems  before 
finishing  with  calibration  of  the  waveguide  system.  Chapter  3  introduces  the 
measurement  setup  and  calibration  scheme,  develops  a  method  for  processing  the  raw 
data  retrieved  from  the  NWA,  and  discusses  how  the  processed  data  is  utilized  by  the 
new  technique  for  extracting  complex  permittivity.  Chapter  4  covers  the  uncertainty 
associated  with  the  measurement  system,  presents  data  gathered  from  measuring 
Plexiglas  samples  of  different  thicknesses,  uses  the  processed  data  to  extract  relative 
complex  permittivity  and  compares  the  results  against  known  standards.  Chapter  5 
contains  conclusions  as  well  as  recommendations  for  future  work. 


1-4 


2.  Background 


Since  the  probe  used  in  this  research  involves  a  parallel-plate  structure,  a  field 
analysis  of  this  guide  is  required  to  understand  the  wave  propagation  behavior  which  is 
necessary  for  the  material  extraction  technique  introduced  in  this  work.  First,  Maxwell’s 
Equations  are  reviewed,  and  a  specialized  formulation  for  this  research  is  given.  Then  a 
Green’s  function,  vector  potential  technique  is  briefly  reviewed.  Next  the  practical  aspect 
of  probe  calibration  is  discussed  due  to  its  importance  in  obtaining  accurate/reliable 
experimental  measurements  for  comparison  to  the  Green’s  function  based  theoretical 
parameter  extraction  model.  Lastly,  the  key  points  of  this  chapter  are  summarized. 

2.1  Maxwell’s  Equations  and  Constitutive  Parameters 

The  derivation  for  this  section  is  drawn  from  Balanis  [4]  and  Harrington  [12]. 
Maxwell’s  equations  in  time  harmonic  form  with  an  assumed  and  suppressed  e’""  time 
dependence  are  given  by 


V  x  E  =  -Mi  -  j'eoB 

(2-1) 

VxH  =  T  +  JC  +  j 

(2-2) 

V*D  =  qev 

(2-3) 

v  •  B  =  qmv 

(2-4) 
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where  E  represents  the  electric  field  intensity;  H  is  the  magnetic  field  intensity;  D  is 


the  electric  flux  density;  B  is  the  magnetic  flux  density;  J  is  the  electric  current  density, 

and  M  is  the  magnetic  current  density.  In  both  cases,  the  subscript  “i”  represents  the 
impressed  or  source  current,  while  with  electric  currents,  the  subscript  “c”  represents  the 
conduction  current  density.  Finally,  qev  is  the  electric  charge  density  while  qmv  is  the 
magnetic  charge  density. 

These  equations  become  well  posed  upon  utilizing  the  constitutive  relations  for 
simple  ( i.e.,  linear  ,  homogeneous,  isotropic)  media,  namely 

D(r,o)  =  s(co)E(r,co)  (2-5) 

B  (r ,  co )  =  //  (®)H(r ,  co) 

Jc  (r ,  co)  =  cr(co)  E  (r ,  co). 

The  previously  determined  scope  of  a  low  loss (e"  «  s'),  non-magnetic  (//(co)  ju0), 

low  dispersive  (s(co)  — >  e)  material  leads  to  the  following  specialized  fonn  of  Maxwell’s 
Equations  which  will  be  utilized  in  this  research 

VxH  =  Ji+crE  +  jcos  E  (2-6) 

=  Ji  +  jco 


s 


c 


.a 

s-J  — 

V  oo) 


—  £q(£ 


(2-7) 


2.2  Green’s  Function  and  Vector  Potentials  in  a  Parallel-Plate  Region 

The  rigorous  formulation  of  the  transverse  fields  inside  the  parallel  plate  region  is 
found  by  replacing  the  waveguide  aperture  with  a  magnetic  current  in  accordance  with 
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Love’s  equivalence  principle  [11].  This  transverse  magnetic  field  can  be  found  by 
utilizing  the  following  relation  [4],  [8] 

HfP  =  — - — {k2  +  V;V  *)f  .  (2-8) 

jcos/u 

where  V,  is  the  transverse  gradient  operator, 


„  d  „  5 
V.  =  x —  +  y — 
dx  dy , 


and  the  electric  vector  potential  F  is  given  as 

b  a 

F  =  ||  G(x,  y,  z  I  x',  y',0)«  <?M(x',  y')dx'dy' 

0  0 


(2-9) 


(2-10) 


where  b  and  a  are  the  aperture  dimensions,  and  G  is  the  rectangular  form  of  the  dyadic, 


parallel-plate  Green’s  function.  G  is  derived  by  Hanson  and  Yakovlev  [5],  [8]  as 


G  =  xGtx  +  yGty  +  zGtz 


-*  uu 

G.„  =  7 -IT  j  f g.AZmz I  z'y^e'^d^n 

(2  n) 

cosh  p  {cl- 1  z  -  z'  l)  ±  cosh  p  {cl- 1  z  +  z'  l) 


S  t  n 


2  p  sinh  pel 


(2-11) 


(2-12) 


(2-13) 


where  d  is  the  sample  material  thickness;  d,  and  //  are  the  Fourier  pairs  for  x  and  y 
respectively;  the  subscript  “f”  and  “n”  identify  the  transverse  and  nonnal  components  of 
the  Green’s  function  respectively;  and  the  z-directed  wave  number  for  the  rectangular 


form  is  p  =  +r f  - k 2  [8]. 
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However,  in  order  to  gain  full  insight  into  the  wave  behavior  in  the  parallel-plate 
region,  the  polar  form  must  be  examined.  Substituting  x  —  x'  =  R cos  6 ,  y  —  y '  =  R  sin  0 , 
i;  =  /l  cos  ([> .  and  q  =  Asin  (/>  into  (2-12)  yields 


-.oo  n 

G  t,n  \eJ 

(2 n)  J0 


JAR  cos(^- 


(2-14) 


where  A  =  ^2  +  q2  and  the  distance  from  the  source  to  the  observation  point  is 

R  =  \p  —  p'\ .  Furthermore,  the  z-directed  wave  number  is  now  defined  as  p  =  yl A2  -k2 
.  Closer  examination  of  the  integral  with  respect  to  (j)  reveals  it  to  have  a  well  known 
solution  of  2tiJ q{AR).  Using  that  solution  with  the  Bessel-Hankel  Function  relation, 

Jo  (z)  =  ~  [//o"  t)  +  Ho2>  (z)]>  and  the  analytic  continuation  formula,  H{2>  (z)  =  -H'2’  (-  z), 
(2-14)  can  now  be  defined  as 


^  oo 

G„  =— J  Xg,p-z\z')H«'(XR)dX 

—00 


(2-15) 


Using  the  complex  plane  analysis  found  in  [8],  the  polar  form  for  (2-15)  is  found  to  be 
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where  l  is  the  parallel-plate  mode  number;  A,  =  ■yjk2  -( iTv/d )2  ;  and  Sl  0  =  1  for  /  =  Oand 
SlQ  -  0  for  /  ^  0  .  Realizing  that  /  =  0  is  both  non-trivial  and  the  dominant  mode, 
several  key  insights  into  the  wave  behavior  inside  the  parallel-plate  region  can  be  made. 
First,  in  the  0th  order  dominate  mode,  the  propagation  constant  is  A0  =  k  =  (o^jpi()  . 


Secondly,  since  the  second  order  Hankel  function  can  be  approximated  in  the  far-field  as 
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-JAR 

the  radial  field  behavior  in  the  far- field  of  the  parallel  plate  region  is  proportional  to  e~jAR 
[8].  Understanding  the  far- field  behavior  is  necessary  to  the  development  of  this  method 
since  an  accurate  characterization  of  the  edge  reflection,  which  is  considered  to  occur  in 
the  far-field,  is  necessary  for  its  use  in  parameter  extraction. 

2.3  Calibration  of  the  Network  Analyzer 

In  order  to  take  accurate  measurements  using  a  parallel  plate  system,  rather 
complicated  calibration  is  often  necessary  [3].  Calibration  complexity  is  increased  by  the 
use  of  a  rectangular  waveguide  probe  which  is  non-standard  for  the  NWA.  Therefore,  in 
order  to  properly  design  a  calibration  method  for  the  rectangular  waveguide  probe,  it  is 
necessary  to  understand  how  the  more  common  calibrations  are  accomplished. 

In  order  to  fully  correct  for  all  of  the  errors  that  the  NWA  may  detect  within  a 
system  it  is  necessary  to  perform  what  is  known  as  a  full  two-port  calibration.  Typically, 
this  is  accomplished  via  one  of  two  methods  utilizing  different  calibrations  standards. 

The  first  is  the  Short-Open-Load-Thru  (SOLT)  method,  and  the  second  is  the  Thru- 
Reflect-Line  (TRL)  technique.  The  latter  (which  is  employed  in  this  research)  is 
commonly  used  in  noncoaxial  systems  due  to  the  ease  of  fabrication  for  the  standards  [6]. 

The  mathematical  model  that  relates  the  S-parameters  of  the  device  under  test 
(5/j,  S2\,  5*22 ,  Sy 2)  to  the  S-parameters  measured  by  the  NWA  detectors  (5™,  5™, 

S22  ,  .S',"")  is  given  by  [6] 
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A  careful  examination  of  these  equations  reveals  that  there  are  12  unknowns  (5,', ,  S2A2 , 
S  A2 ■  $n ,  S2*2 ,  5n  ,  5,c; ,  52c2 . SAS  A  ,  ,  5f252 j ,  and  S*S£ ) .  In  the  TRL  technique 
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and  Line 
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standard  (in  addition  to  other  assumptions)  are  utilized  to  solve  for  the  error  terms  in  the 
calibration  model.  Once  these  error  terms  are  determined,  the  desired  S-parameters  of 
the  sample  measurement  can  be  accurately  computed  using  the  NWA  processor. 

2.4  Summary 

In  this  chapter,  a  brief  review  of  EM  theory  was  provided.  Additionally,  the  wave 
behavior  inside  a  parallel  plate  system  was  discussed.  This  characterization  provides  the 
foundation  of  the  mathematical  assumptions  made  in  Chapter  3.  Most  of  the 
characterization  was  limited  to  the  very  specific  example  of  parallel  plate  regions  as  this 
is  a  vital  part  of  the  measurement  system  used  in  this  thesis.  More  detailed  information 
on  wave  behavior  and  Green’s  functions  can  be  found  in  the  sources  cited.  Finally,  a 
basic  TRL  calibration  method  was  discussed  in  order  to  provide  the  background  of  the 
more  specialized  TRL  calibration  that  will  be  developed  in  Chapter  3. 
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3.  Methodology 


Based  on  the  information  presented  in  previous  chapter,  a  method  for  extracting 
permittivity  from  a  non-magnetic,  non-dispersive,  low  loss  dielectric  using  a  waveguide 
probe  is  detailed.  The  process  includes  setting  up  and  calibrating  the  measurement 
system,  determining  the  frequency  range,  processing  the  data,  and  concludes  with  the 
mathematical  extraction  of  the  permittivity  from  the  measured  and  processed  sample  S 
parameters. 

3.1  Measurement  System  Setup 

The  measurement  system  consists  of  several  parts,  as  shown  in  Figure  3-1.  The 
system  includes  a  NWA,  cables,  Ku-band  waveguides,  and  flange  plates  to  be  attached  to 
the  end  of  the  waveguides.  Note,  although  a  single  probe  is  used  for  material  property 
extraction,  a  full  2-port  TRL  calibration  is  performed. 

All  flange  plates  contain  a  Ku-band  sized  aperture  that  allows  for  fields  from  the 

waveguide  to  pass  through  into  the  parallel-plate  sample  region.  The  flange  plates  are 

split  into  two  separate  geometries  (circular  and  square)  and  have  two  different 

dimensions  per  geometry.  For  the  circular  flange  plates,  the  diameters  are  6  inches 

(0.1524m)  and  4  inches  (0.1016m).  For  the  square  flange  plates,  the  larger  ones  have  a 

length  and  width  of  6  inches  (0.1524m)  while  the  smaller  ones  have  a  length  and  width  of 

4  inches  (0.1016m).  The  square  flange  plates  mimic  geometries  used  in  previous  parallel 

plate  waveguide  probe  work  [9].  However,  the  circular  flange  plates  were  developed  in 

an  attempt  to  increase  the  amplitude  of  the  response  from  the  edge  reflection.  For  each 
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Figure  3-1:  Measurement  system  components:  (a)  Network  Analyzer,  (b)  Ku- 
band  waveguide,  (c)  4”  length/width  square  plates,  (d)  6”  length/width  square 
plates,  (e)  cable  to  connect  waveguide  to  NWA,  (f)  4”  diameter  circular  plates, 
and  (g)  6”  diameter  circular  plates. 


plate  size  and  geometry,  there  are  two  plates  in  order  to  facilitate  a  two  port  TRL 
calibration.  Drawings  for  the  plates  are  located  in  Appendix  B. 

In  order  to  setup  the  measurement  system,  a  flange  plate  is  attached  to  the  end  of 
the  waveguide  so  that  the  apertures  line  up  (using  precision  alignment  pins),  and  the 
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Figure  3-2:  Fully  connected  waveguide  probe. 

cable  is  used  to  connect  the  probe  to  the  NWA.  Figure  3-2  shows  a  connected 
waveguide  probe. 

3.2  Calibration  Details 

As  discussed  in  the  previous  section,  there  are  two  plates  of  each  size  and 
geometry  in  order  to  facilitate  a  two  port  calibration.  The  two  port  calibration,  while 
more  intensive  than  no  calibration  or  a  one  port  calibration,  carries  many  benefits.  First, 
a  two  port  calibration  identifies  and  removes  errors  from  the  measurement  system. 
Additionally,  since  the  system  is  waveguide  based,  the  TRL  method  can  be  used.  This 
carries  the  benefit  of  easy-to-manufacture  standards  [6] .  Finally,  and  most  importantly, 
the  two  port  calibration  allows  for  the  reference  plane  to  be  set  at  the  point  where  the 
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Figure  3-3:  Calibration  setups  for  the  (a)  thru  standard,  (b)  line  standard, 
and  (c)  reflect  standard. 

plate  aperture  meets  the  sample  material.  This  means  that  all  temporal  reflections 
measured  within  the  probe  system  will  be  referenced  to  time  t0  =  0  . 

However,  because  of  the  unique  geometry  associated  with  the  parallel  plate 
system,  the  two  port  TRL  calibration  has  to  be  modified  slightly.  The  thru  standard, 
shown  in  Figure  3-3(a),  is  formed  by  joining  the  probes  connected  to  each  port  of  the 
NWA  so  that  there  is  a  single  open  path  from  one  port  to  the  other.  With  the  two 
calibration  planes  of  each  probe  joined  up,  there  will  be  no  reflections  at  the  sample. 
Therefore,  the  S2’,  and  Sjs2  values  will  be  1  since  all  the  energy  transmitted  down  one 
waveguide  will  be  transferred  to  the  other.  This  coincides  with  the  sample  setup  in 
(2-19). 
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Figure  3-4:  Diagram  showing  the  size  difference  parallel  plate  thru  standard 
setup  (denoted  by  solid  lines  with  a  dash  line  separating  the  two  probes)  and  the 
line  standard  (denoted  by  the  dotted  lines). 

The  line  standard  presents  a  tougher  problem  for  calibration  since  parallel  plate 
systems  are  designed  to  be  non-destructive.  However,  the  calibration  measurement  can 
be  accomplished  by  using  the  regular  waveguide  line  standard,  shown  in  Figure  3-3(b). 
Examining  (2-21)  gives  insight  into  how  this  is  possible.  In  a  regular  waveguide 
calibration,  the  line  standard  measurement  is  calculated  by  using  the  distance  between  the 
two  calibration  planes  as  an  offset  in  order  to  calculate  the  field’s  phase  changes.  Using 
this  type  of  standard,  the  distance  is  a  positive  value.  However,  as  is  evident  from  Figure 
3-4,  the  distance  between  the  two  reference  planes  is  no  longer  positive.  Therefore,  if 
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this  negative  distance  can  be  calculated,  then  the  negative  phase  change  can  be  calculated 


using  the  formula 


T  _  —  d  _  ^ plates  d LS  (3-1) 

c  c 

where  dplates  is  the  size  of  two  plates  stacked  together  and  dLS  is  the  size  of  the  regular 

waveguide  line  standard.  The  line  standard  used  for  the  measurement  system  has  a 
thickness  of  6.44mm.  Using  this  value  in  conjunction  with  (3-1),  the  delay  time  for  each 
plate  size  and  geometry  can  be  calculated  as  shown  in  Table  3-1. 


Table  3-1:  Plate  thickness  measurements  and  associated  delay  time 


Meas  1 
(mm) 

Meas  2 
(mm) 

Meas  3 
(mm) 

Meas  4 
(mm) 

Avg  (mm) 

d  (mm) 

M  ps) 

6"  Rect 

13.14 

13.23 

13.13 

13.12 

13.155 

6.715 

-22.3988 

4"  Rect 

13.15 

13.09 

13.15 

13.12 

13.1275 

6.6875 

-22.3071 

6"  Circ 

13.17 

13.12 

13.10 

13.13 

13.13 

6.69 

-22.3154 

4"  Circ 

13.16 

13.14 

13.14 

13.15 

13.1475 

6.7075 

-22.3738 

The  reflect  standard  is  easily  realized  via  placement  of  a  metal  plate  between  the 
probe  flanges  as  shown  in  Figure  3-3(c)  This  coincides  with  the  setup  given  in  (2-20). 

As  a  final  note,  when  accomplishing  the  TRL  calibration,  in  order  to  minimize  the 
flange  plate  reattachment,  the  calibration  is  accomplished  by  doing  the  line  standard  first, 
then  the  thru  and  reflect  standards. 


3.3  Frequency  Range 

Prior  to  calibrating  the  NWA,  the  frequency  range  for  the  measurement  must  be 
set.  The  waveguide  probe  contains  two  structures  that  present  limits  to  the  frequency 
range.  The  first  is  the  waveguide  itself.  A  full  discussion  of  the  frequency  range  that  the 
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Ku-band  waveguides  used  in  this  research  can  cover  is  presented  in  Appendix  C.  The 
second  structure  is  the  parallel-plate  region  created  by  the  flange  plate  and  the  PEC 
backed  sample.  The  cutoff  frequency,  fc ,  for  parallel  plate  regions  like  these  is  given  as 


where  n  is  the  mode  number  and  a  is  thickness  of  the  sample  [4].  For  the  0th  order  mode, 
which  as  discussed  in  Chapter  2  is  dominant,  fc  =  0Hz .  However,  depending  on 
complex  permittivity,  sample  thickness,  and  frequency,  the  1st  order  mode  can  be  excited 
in  the  parallel-plate  region.  Once  “higher  order”  modes  are  excited,  the  conclusions  from 
Chapter  2  about  wave  behavior  in  the  parallel-plate  region  are  no  longer  valid.  Since 
sample  type  and  thickness  may  be  fixed,  either  the  theory  should  be  modified  or  the 
upper  limit  of  the  frequency  range  for  the  system  should  be  reduced. 

Furthermore,  (3-2)  requires  the  relative  complex  permittivity  in  order  to 
determine  cutoff  frequency.  This  presents  a  problem  since  sr  is  what  is  supposed  to  be 
determined  by  the  method.  As  a  solution  to  this  problem,  a  general  estimate  for  the 
cutoff  frequency  can  be  obtained  heuristically  by  observing  when  the  frequency-domain 
Su(co)  measurement  begins  to  oscillate  heavily  and  no  longer  behaves  in  a  manner 
similar  to  the  lower  frequency  part  of  the  measurement,  that  point  can  be  estimated  to  be 
the  cutoff  frequency.  An  example  of  this  is  shown  in  Figure  3-5. 

3.4  Data  Processing 

Upon  completion  of  system  calibration,  conductor-backed  sample  measurements 
can  be  performed.  This  is  accomplished  by  placing  a  PEC  backed  sample  material  under 
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Figure  3-5:  Frequency-domain  data  from  a  waveguide  probe  that  has  the  1st 
order  mode  being  excited  in  the  parallel-plate  region 

the  probe  in  a  setup  similar  to  Figure  3-2  and  performing  a  frequency  swept  NWA  Su(eo) 
measurement.  The  temporal  response  can  be  subsequently  obtained  via  inverse  Fourier 
transformation  (i.e.,  IFFT)  of  the  spectral  domain  response,  Sn(co)  .  However,  the 
response  from  the  edge  of  the  plate  may  not  be  distinguishable  from  the  side  lobes  of  the 
main  reflection.  An  example  of  this  is  shown  in  Figure  3-6.  Therefore,  additional  signal 
processing  is  necessary  in  order  to  extract  the  edge  reflection. 

In  order  to  properly  process  the  data,  the  transform  processes  used  by  the  NWA 
must  be  fully  understood  in  order  to  accurately  extract  the  edge  reflection  from  the  data. 
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Figure  3-6:  Raw  time  domain  data  of  a  4”  diameter  circular  plate 
3.4.1  Fourier  Analysis 

A  common  operation  performed  in  signal  processing  is  the  Fourier  transform. 
Due  to  the  discrete  sampling  of  the  measured  NWA  data,  transforms  are  performed  via 
the  discrete  Fourier  transform  (DFT).  The  DFT  transform  pair  is  given  as 


N- 1 

x(n)  ='s^jX\k\e  N 


k=  0 


.2  n , 


r  n  1  "4  ,  ,  -j—kn 

xik]  =  T;jlx(n)e  N 

N  n= 0 


(3-3) 


(3-4) 
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where  N  is  the  number  of  samples  for  both  the  frequency  domain  data,  X\k\  and  the 
time  domain  data,  x{n ) .  Due  to  the  inherent  properties  (i.e.,  periodic  nature)  of  the 
DFT,  it  is  of  critical  importance  to  understand  the  effects  that  finite  spectral  bandwidth 
( B  W)  and  frequency  step  size  (  df)  have  on  temporal  resolution  (Ts)  and  alias-free  time 
span  (TAf).  It  is  also  important  to  understand  the  effect  of  zero-padding  data  (which  is 
employed  to  achieve  2N  data  points  required  by  the  computationally  effective  FFT 
algorithm).  Using  the  information  above,  the  temporal  resolution,  7) ,  is  approximately 
determined  using  the  frequency  step  size,  df  ,  via 

7>  — =  —  (3-5 

Ndf  BW 

The  alias-free  range  for  the  time-domain  data  can  be  determined  via  TAF  =NTS  [14]. 

Using  the  relations  given  above,  frequency  data  collected  from  the  NWA,  shown 
in  Figure  3-7 (a),  can  be  transformed  into  time-domain  data,  shown  in  red  in  Figure 
3-7 (b),  by  using  the  transform  relations  of  (3-3)  and  (3-4).  When  comparing  the 
transformed  data  with  data  gathered  from  the  NWA,  plotted  in  blue  in  Figure  3-7(b),  it 
becomes  apparent  that  the  two  graphs  do  not  match.  However,  as  previously  stated,  the 
data  gathered  from  the  NWA  is  not  directly  measured  but  the  product  of  a  Fourier 
analysis  performed  by  the  NWA  on  the  frequency-domain  data. 
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Figure  3-7:  (a)  Measured  frequency  domain  data  and  (b)  data  transformed  via  the 
NWA  (NWA  data)  and  data  transformed  using  the  DFT  relations  of  (3-3)  and 
(3-4)  (FT  data)  from  a  probe  using  4  inch  circular  plate. 

This  means  that  some  other  process  is  at  work  in  the  NWA’s  Fourier  analysis. 

The  most  likely  reason  for  the  apparent  increase  in  resolution  shown  in  the  NWA  data  in 
Figure  3-7(b)  is  a  process  known  as  zero  padding.  The  NWA  is  setup  to  save  however 
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many  points  of  data  the  user  specifies.  This  is  true  for  both  frequency  data  and  time  data. 
Additionally,  the  NWA  will  also  save  the  number  of  points  for  the  time  range  specified. 

As  an  example,  if  the  user  specifies  that  the  NWA  should  have  401  points  of  data 
and  sets  the  range  for  the  time  data  to  cover  -2  nano-seconds  (ns)  to  2  ns,  then  the  NWA 
will  save  401  data  points  within  that  time  range.  However,  this  presents  a  problem. 

Since  the  NWA  does  not  actually  measure  time  domain  data,  the  data  displayed  has  to  be 
calculated  using  some  variation  on  the  DFT.  This  means  that  the  relation  set  up  in  (3-5) 
applies.  Therefore,  the  existing  number  of  samples  has  to  fit  into  the  time  range 
determined  by  that  relation.  If  this  range  exceeds  the  range  set  by  the  user  in  the 
example,  then,  in  order  for  the  NWA  to  still  be  able  to  save  the  user  defined  data  points, 
it  needs  to  use  zero  padding  to  increase  the  number  of  samples  so  that  the  saved  data  has 
the  user  defined  number  of  points.  This  principle  is  shown  in  Figure  3-8  where  the 
transformed  data  is  zero  padded  until  it  matches  the  data  from  the  NWA.  The  final  graph 
shows  that  when  1601  data  points  are  used  for  the  zero  pad,  the  computed  data  matches 
the  data  saved  from  the  NWA.  The  large  number  of  data  points  required  for  the  zero-pad 
is  necessary  since  the  unpadded  transformed  data  only  had  34  data  points  covering  the 
region  from  -2ns  to  2ns.  After  application  of  the  pad,  the  new  data  set  now  has  401 
points  covering  the  region  from  -2ns  to  2ns. 
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Figure  3-8:  Comparison  time-domain  data  processed  by  the  NWA  (shown  in 
blue)  with  data  processed  using  the  DFT  algorithm  (shown  in  red)  with  a  zero 
padding  of  (a)  0,  (b)  101  and  (c)  1601  points  for  a  probe  with  a  4  inch  circular 
plate. 

3.4.2  Kaiser  Windowing 

After  understanding  the  signal  processing  accomplished  by  the  NWA,  additional 
signal  processing  is  necessary  for  extracting  the  edge  response.  The  standard  processing 
tool  available  on  the  NWA  is  a  Kaiser  window.  The  Kaiser  window  was  chosen  in  this 
research  to  represent  the  effect  of  windowing  functions.  This  decision  was  made  since 
the  research  is  designed  to  be  a  proof  of  concept,  the  application  of  the  window  will  give 
a  general  idea  of  how  similar  windows  affect  the  data.  Additionally,  since  the  Kaiser 
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window  is  a  built  in  function  of  the  NWA,  using  it  maintains  the  simplicity  of  the  method 


as  no  new  processing  algorithms  are  needed. 
The  Kaiser  window  is  defined  as 


w(t)  =  <! 


Wl-(2  tf 


hiP) 


I  I  1 

t  <  - 
1  1  2 


elsewhere 


(3-6) 


in  the  time  domain  and 


,f(/)  sin(A72-l2)  (3-7) 

in  the  frequency  domain.  In  both  domains,  ft  covers  a  range  from  0  to  10  [10].  In  the 
Kaiser  window,  70  is  a  0th  order  modified  Bessel  function. 
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Figure  3-9:  Time-domain  data  set  with  (a)  no  Kaiser  windowing  and  with  the 
following  Kaiser  windows  applied  (b)  fi=l,  (c)  /J-3,  and  (d)  p=6. 

The  NWA  can  apply  the  Kaiser  window  on  using  an  adjustable  /?  value  in  order 
to  suppress  the  sidelobes  of  the  main  response  so  that  smaller  responses,  like  the  edge 
reflection,  become  visible.  An  example  NWA  data  set  with  no  zero  padding  and  varying 
[j  values  is  shown  in  Figure  3-9.  While  waveform  representing  the  edge  reflection  is 
easily  identifiable,  waveform  spreading  is  noticed.  Furthermore,  the  peak  of  the  edge 
reflection  moves  with  the  application  of  the  window.  This  offset  is  somewhat  mitigated 
through  the  use  of  zero-padding  applied  by  the  NWA.  Therefore,  when  examining  data 
with  a  Kaiser  window  applied,  the  results  will  be  zero-padded. 
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3.4.3  Signal  Subtraction 

Another  method  for  extracting  the  edge  reflection  takes  advantage  of  the  linearity 
of  the  system  by  using  the  principle  of  superposition  to  isolate  the  desired  data.  When 
superposition  applies,  it  dictates  that  the  signal,  in  either  the  frequency  domain  or  the 
time  domain  can  be  broken  into  smaller  parts  that  can  be  added  together  to  form  the 
original  signal  [14].  Therefore,  by  utilizing  physical  intuition  about  the  system  setup  and 
how  it  will  respond  when  excited,  certain  approximations  can  be  employed  to  aid  in 
indentifying  the  edge  reflection  signal. 

The  first  step  is  to  analyze  the  system  and  determine  what  is  physically  happening 
when  it  is  excited.  Figure  3- 10(a)  depicts  that  the  EM  waves  reflect  off  of  the  sample  and 
the  PEC  backing  while  other  EM  waves  move  outward  away  from  the  center.  These 
outgoing  waves  will  reflect  off  of  the  edge  of  the  plate,  yet  not  all  of  the  waves  will 
reflect  back  toward  the  center.  Some  waves  will  continue  to  head  outward  into  the 
sample  material.  Furthermore,  the  waves  will  undergo  some  attenuation  as  they  travel 
outward  toward  the  edge  and  back.  This  is  shown  in  further  detail  in  Figure  3-11. 

Taking  this  information  into  account,  it  is  relatively  easy  to  deduce  that  the  strongest 
response  (the  0th  order  reflection)  will  be  from  the  initial  reflection  off  of  the  PEC 
backing  of  the  sample  while  other  responses,  like  the  edge  reflection  (the  1st  order 
reflection),  will  not  be  as  strong  when  they  return  to  the  probe. 
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PEC  Surface 

(a) 


(b) 


Figure  3-10:  (a)  Drawing  showing  simplified  wave  behavior  in  the  measurement 


system  and  (b)  measured  frequency  data  from  a  probe  using  a  4  inch  plate. 
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To  VNA 


Figure  3-11:  Diagram  depicting  the  material  reflection  Su,o,  the  first-order 
edge  reflection  Su,u  and  the  second-order  edge  reflection  Sn,2- 

With  this  physical  insight  in  mind,  Figure  3- 10(b)  can  be  examined. 

Superposition  says  that  the  individual  frequency  responses  from  each  scattering  point  will 
add  together  into  the  final  signal  in  the  following  manner 

*^ii  (®)  —  *^1 1 jo  (®)  ^11,1  (®)  (®)  •  •  (3-8) 

where  St  M  (a)  is  the  frequency-domain  measured  data  and  SUfl  (co)  is  the  nth  order 

reflection.  Therefore,  if  there  is  one  overly  dominant  response,  like  the  0th  order 
reflection,  it  will  drive  the  basic  shape  of  the  response  curve  with  the  smaller  reflections 
being  represented  as  minor  contributions.  Based  on  the  conclusion  that  the  0th  and  1st 
order  reflections  are  much  stronger  than  the  higher  order  reflections,  (3-8)  can  be 
approximated  as 
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snM  O)* 

5np  (co)  +  Sn  i  (co) 


(3-9) 


Furthermore,  since  Sn  f)  (co)  »  Sj  Ll  (co)  ,  if  the  basic  shape  of  Sj  "  ( co )  can  be  recreated 

via  some  kind  of  low-order  polynomial  fit  without  accounting  for  the  minor 
contributions,  the  0th  order  reflection  can  be  approximated.  The  0th  order  reflection,  once 

approximated,  can  be  subtracted  from  SnM  (co) .  Based  on  (3-9),  this  should,  in  effect, 
leave  the  only  the  1st  order  response  left. 

Figure  3-12  illustrates  how  this  principle  works.  In  Figure  3-12(a),  a  low-order 
polynomial  fit  is  applied  to  measured  frequency-domain  data.  Based  upon  physical 
intuition  and  superposition,  this  fit  should  be  an  approximation  of  the  0th  order  reflection. 
In  Figure  3-12(b),  a  transform  is  applied  to  the  fit,  and  when  compared  to  the  NWA  time 
domain  data,  the  fit  serves  as  a  close  approximation  of  the  0th  order  reflection.  It  should 
be  noted  that  the  transform  conducted  on  the  fit  did  not  use  zero  padding.  If  the  fit  in 
Figure  3- 12(b)  is  compared  with  the  Fourier  analysis  done  on  the  frequency  data  in 
Figure  3-8(a),  the  approximation  is  much  closer.  Finally,  after  approximating  the  0th 
order  reflection,  utilizing  superposition,  the  reflection  is  subtracted  from  the  time  domain 
data.  This  is  shown  in  Figure  3- 12(c).  Once  the  0th  order  reflection  is  removed,  the  only 
major  response  left  is  the  1st  order  reflection  which  is  the  edge  reflection.  Thus  the 
desired  delay  time  is  now  easily  discernable. 
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Figure  3-12:  Process  of  performing  a  dominant  signal  subtraction:  (a)  a  polynomial 
fit  is  applied  to  the  frequency  domain  data,  (b)  the  fit  is  transformed  into  an 
approximation  of  the  0th  order  reflection,  and  (c)  the  response  is  subtracted  from  the 
time-domain  data  leaving  the  edge  reflection. 
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3.5  Extracting  values 

After  isolating  the  edge  response  so  that  the  time  of  its  occurrence  and  amplitude 
can  be  determined,  the  permittivity  can  be  calculated.  However,  in  order  to  develop  a 
simple  extraction  technique,  some  assumptions  are  made  based  upon  the  formulas 
reviewed  in  Chapter  2. 

Based  upon  the  conclusions  reached  via  (2-16),  the  propagation  constant,  k  ,  for 
the  waves  inside  the  parallel  plate  structure  is  defined  as  k  =  (o^qx0  =  ft  -  ja  where  [3 
defines  the  propagation  factor  and  a  defines  the  attenuation  factor.  The  phase  velocity, 
vp  ,  for  the  wave  can  be  calculated  via  the  formulation 


CO  _  CO 
p  ~  Re(k)  ~  ~p 


(3-10) 


where  co  =  2nf  is  the  angular  frequency.  Based  on  the  information  given  in  (3-10)  and 
(2-7),  (3  and  a  can  be  defined  as 

P  =  ReW  =  Re(®V^0)  =  Rc(W £oMo  a  )  =  Rc(k»  ^ )  =  R etk^e'-je")  (34 1) 


a  =  -Im(/:)  =  -Im  (kQ^e'-je”) . 


(3-12) 


Since  velocity  is  the  ratio  of  distance  over  time,  and  given  that  the  phase  velocity  of  the 
waves  is  defined  in  (3-10),  the  time  it  takes  for  an  EM  wave  to  travel  any  distance  in  the 
parallel  plate  system  can  be  defined  as 


t  = 


v 


p 


(3-13) 
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Thus,  two-way  travel  time  to  the  edge  and  back,  t2w,  which  is  import  for  permittivity 
extraction,  is  easily  computed  via  the  formula 

t,  (3-14) 

2  w 

vp 

where  d2w  is  the  distance  from  the  center  of  the  flange  plate  aperture  to  the  flange  plate 
edge  and  back. 

Another  relation  required  for  permittivity  extraction  can  be  obtained  by 
examining  the  attenuation  that  occurs.  Using  the  previous  assumption  about  wave 
behavior  at  the  edges  of  the  parallel  plate  structure,  the  response  from  the  edge  can  be 
modeled  as 


e 


-adt 


e 


- m 


(3-15) 


where  R  is  a  value  representing  the  complex  magnitude  of  the  reflected  energy  from  the 
edge.  Since  the  amount  of  energy  being  lost  at  the  edges  is  not  known,  R  cannot  be  easily 
found.  However,  if  the  magnitude  of  the  edge  response  from  two  different  size  plates 
were  divided  by  each  other,  then  R  would  cancel  along  with  yfk  ,  leaving  only  a 
unknown  and  therefore  solvable.  Putting  the  two  measurements  together  yields  the  ratio 
A,  which  is  defined  as 


A  - 


(3-16) 


This  is  second  equation  necessary  for  finding  complex  permittivity. 
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3.5.1  Simplified  Extraction 

Now  that  the  necessary  equations  have  been  defined,  in  order  to  compute  complex 
permittivity  in  a  simple  manner,  one  last  assumption  is  made.  This  assumption  uses  the 

binomial  expansion  to  find  a  closed  form  solution  since  s/1  will  not  provide  one. 

(l  +  jc)"  «  1  +  h.v...  x  «  1.  (3-17) 


Applying  this  approximation  to  the  relative  permittivity  in  the  propagation  constant,  k, 
yields 


./  ./  ■./(  C-"V2  ,  .  , /(  p" 

e/1  =  (e'-je'  Y  =  YY2  1  -  J  ~  ~  YY  1  -  j  — 

\  £  )  \  2  e  j 


\Yi 


(3-18) 


Therefore, 


Im(s Y )~~  S 


iYs' 
Re(s/2) «  Ys' 


(3-19) 


which  implies  that 


a  =  -k0  Im(f  /2 )  «  -kr 


f  e"  A 


k0e " 


iYs'  J  iYs' . 


P  =  kQ  R e(sY2 )  «  kQ  Vf7 


(3-20) 


Substituting  (3-20)  into  (3-14)  and  (3-16)  gives 

.  _  d2w  _  d2w  _  Pd2w  ~  kQ  Ys  d2w  (3-21) 

Vp  G)/p  0)  0) 

ln(A)  -x~k"£" 

d2w  -  d2w  ~  2YY 

which  can  be  manipulated  to  solve  for  the  relative  complex  permittivity  with  the  solution 
for  s'  given  as 
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£ 


f 


(3-22) 


cot 


2  w 


\^Qd2w  J 


and  e"  given  as 


2V^ln(A) 

/  I  1 2 w  1 2w  I  * 

h\d2  -d\  ) 
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3.6  Summary 

This  chapter  has  laid  out  a  process  for  setting  up  and  calibrating  a  rectangular 
waveguide  probe.  Furthermore,  frequency  range  considerations  were  presented  in  order 
to  prevent  the  generation  of  “higher  order”  modes  in  the  parallel  plate  region  of  the  probe 
setup.  The  effect  of  this  is  seen  when  applied  to  Plexiglas  samples  in  Chapter  4.  Two 
signal  processing  methods  were  also  discussed,  and  both  are  used  in  the  following 
chapter.  Finally,  the  extraction  technique  was  developed  based  upon  the  formulation 
from  Chapter  2.  The  results  from  using  this  technique  are  presented  in  Chapter  4. 
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4.  Analysis  and  Results 


This  chapter  takes  the  method  laid  out  in  Chapter  3  and  presents  the  data  gathered 
from  using  this  method  on  two  Plexiglas  samples.  First,  the  uncertainty  associated  with 
the  measurements  is  analyzed  so  that  confidence  intervals  for  the  final  extraction  of 
permittivity  can  be  determined.  Then  the  results  of  the  two  signal  processing  methods 
discussed  in  Chapter  3  are  presented.  Finally,  the  permittivity  extraction  technique  is 
applied  to  the  data,  and  results  are  presented. 

4.1  Error  Analysis 

The  material  characterization  technique  developed  in  this  research  will  have  some 
degree  of  error  associated  it.  Understanding  this  error  is  important  in  being  able  to  assess 
the  method’s  accuracy  and  effectiveness. 

Before  the  uncertainty  of  the  relative  complex  permittivity  can  be  computed,  the 
sources  of  error  must  be  identified.  The  formulas  in  (3-21)  provide  the  basis  for 
identifying  these  error  terms.  Of  these  values,  four  are  variable  and  must  be  measured. 
One  of  these,  the  variable  representing  distance  from  the  center  of  the  aperture  to  the 
edge  and  back  ,  d2w ,  is  determined  by  using  calipers  to  measure  the  flange  plates.  The 
other  three  (angular  frequency,  co,  two  way  time,  t2w,  and  the  ratio  of  response 
amplitudes,  A)  are  measured  by  the  NWA,  and  two  of  these  (t2w  and  A)  are  derived 
through  the  Fourier  analysis  accomplished  by  the  NWA. 
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The  NWA  utilizes  averaging  and  IF  bandwidth  to  reduce  measurement  noise 
level.  Since  the  value  A  is  a  ratio,  and  the  signal  processing  will  apply  similar 
mathematical  mechanisms  on  the  signal,  the  effect  of  the  uncertainty  in  A  is  assumed  to 
be  negligible. 

4.1.1  Time  and  Frequency  as  Variables 

The  manual  [3]  lists  the  minimum  frequency-domain  resolution  for  the  NWA  as 
1Hz.  This  resolution  is  assumed  to  be  the  worst  case  uncertainty  for  a  frequency 
measurement.  Given  that  the  frequency  step  size  is  several  orders  of  magnitude  greater 
than  the  uncertainty,  it  is  assumed  this  uncertainty  is  negligible. 

The  time-domain  uncertainty  consists  of  two  parts.  The  first  is  drawn  from  the 
frequency  uncertainty.  However,  because  the  frequency  uncertainty  is  considered 
negligible,  this  part  of  the  time-domain  uncertainty  is  also  considered  negligible. 

Fourier  analysis  provides  the  second  form  of  time-domain  uncertainty  which  has  a 
larger  effect.  Equation  (3-5)  shows  that  the  time-domain  resolution  is  dependent  on  the 
bandwidth.  Therefore,  it  is  possible  that  the  actual  edge  reflection  lies  within  the 
irresolvable  region.  However,  it  is  assumed  that  this  edge  reflection  peak  is  closer  to  the 
point  that  is  identified  as  the  edge  reflection  than  any  other  point.  This  means  that  the 
actual  response,  x,  will  lie  somewhere  within  an  area  of  -Tj2<x<Tj2  where  Ts  is  the 
temporal  resolution.  This  area  is  centered  on  the  point  that  is  identified  as  the  edge 
reflection.  Furthermore,  since  the  actual  edge  reflection  is  not  dependent  on  the  Fourier 
analysis  but  on  the  material  under  test,  it  could  theoretically  occur  anywhere  within  area 
defined  above  with  equal  probability.  Therefore,  the  probability  density  function  (PDF) 
of  the  time  of  occurrence  for  the  edge  reflection  can  be  defined  as  a  unit  step  function  [1], 
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Figure  4-1:  The  histograms  and  associated  Gaussian  curves  for  the  measurement 
plates 

This  means  that  the  standard  uncertainty,  u,  for  the  time  of  occurrence  of  the  edge 
reflection  can  be  defined  as 

u  =  \  (4-D 

2v3  . 

4.1.2  Distance  as  a  Variable 

The  two  way  distance  for  each  flange  plate  was  measured  using  precision 
calipers.  A  sample  of  20  measurements  per  plate  was  taken  in  order  to  be  able  to  form  a 
good  statistical  picture.  The  plates  were  found  to  exhibit  a  Gaussian,  or  normal,  PDF 
with  a  low  standard  deviation.  The  distributions  of  the  selected  plates  with  an  overlaid 
Gaussian  curve  are  shown  in  Figure  4-1.  The  statistical  information,  shown  in  Table  4-1, 
forms  the  basis  of  the  uncertainty  associated  with  distance  in  the  extraction  formulas. 
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Table  4-1:  Statistical  data  generated  from  20  sample  measurements  of  each  plate 
that  was  used  in  taking  measurements. 


Mean  (mm) 

Standard  Deviation  (mm) 

6”  Circular  Plate 

152.38 

0.0186 

4”  Circular  Plate 

101.59 

0.0145 

6”  Square  Plate 

152.42 

0.0112 

4”  Square  Plate 

101.61 

0.0142 

4.1.3  Overall  Uncertainty 

As  shown  in  the  previous  sections,  both  the  uncertainty  in  the  distance  traveled  by 
the  wave  in  the  parallel  plate  system  and  the  uncertainty  of  the  time  variable  resolution 
have  the  most  significant  effect  on  the  results  of  the  permittivity  extraction.  Therefore, 
the  uncertainty  of  both  must  be  taken  into  account  during  the  permittivity  extraction.  For 
the  real  part  of  the  relative  complex  permittivity,  uncertainty  can  be  computed  as  [13], 
[15] 

8s' (d,  t )  =  d£'(d°  - to ^  Sd  +  d£'(d°  St  (4-2) 

x  ’  dd  dt 


with  d  and  t  being  the  distance  and  time  variables,  the  subscript  0  denotes  the  expected 
values  and  the  S  denotes  the  uncertainty.  The  imaginary  part  of  the  relative  complex 
permittivity  is  defined  similarly,  as 

A  ^ '  \  A^"lA  a  ^ '  \ 

Sd „  (4-3) 


5s 


,  .  ds"(d,  ,d0  ,s')  ds"(d ,  ,d9  ,s') 

”(d],d2,s')= - V  lo  2°  °JScL  + - V  lo  2°  0J 

y  1  2  J  dd{  1  dd2 


ds”\d ,  ,d9  ,s') 

_ V  !()  Z0  U  /  ggf 

ds' 
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where  dt  and  are  the  expected  values  for  dx  and  d2  respectively  and  s'0  is  the 
expected  value  for  s' . 

4.2  Results  of  Data  Processing 

Before  the  data  can  be  used  in  calculations  to  find  the  relative  complex 
permittivity,  it  must  first  be  processed  using  one  of  the  methods  detailed  in  Chapter  3. 
Each  method  carries  with  it  benefits  and  drawbacks.  However,  the  overall  extent  of  the 
drawbacks  remains  unknown  until  the  extraction  is  accomplished  and  can  be  compared 
against  known  values.  As  a  preliminary  measure  of  effectiveness,  each  figure  shown 
below  will  have  a  marker  showing  the  ideal  location  of  the  edge  response  using  the 
industry  standard  value  of  2.6  for  the  real  part  of  relative  complex  permittivity. 

4.2.1  Kaiser  Windowing 

As  a  processing  method,  Kaiser  windowing  presents  an  attractive  option, 
especially  since  the  NWA  has  the  ability  to  perform  the  windowing  on  the  data.  It  is  able 
to  suppress  the  sidelobes  of  the  signals  so  that  weaker  signals,  like  the  edge  reflection, 
can  be  more  readily  identified.  However  this  comes  at  a  cost  of  spreading  the  waveforms 
of  the  responses  as  shown  in  Figure  3-9.  The  figures  below  present  time  domain  data 
that  has  been  taken  from  5.55mm  and  4.39mm  thick  samples  of  Plexiglas.  A  frequency 
ranges  were  10.0  GHz  to  16  .0GHz  and  10.0  GHz  to  18.5  GHz  respectively.  The  thicker 
sample  has  a  reduced  frequency  range  to  avoid  exciting  the  1st  order  mode.  The  Kaiser 
window  was  applied  using  the  NWA  with  a  [1  value  of  6  with  zero  padding.  Data  was 
taken  using  the  4”  and  6”  circular  plates.  As  a  frame  of  reference,  each  figure  below  also 
has  data  using  the  same  setup  without  the  Kaiser  window  plotted  with  it. 


4-5 


Figure  4-2:  Data  from  the  6”  circular  plate  measuring  a  5.55mm  thick  sample 


Figure  4-3:  Data  from  the  4”  circular  plate  measuring  a  5.55mm  thick  sample 
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Figure  4-4:  Data  from  the  6”  circular  plate  measuring  a  4.39mm  thick  sample 


Figure  4-5:  Data  from  the  4”  circular  plate  measuring  a  4.39mm  thick  sample 
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As  is  evident  from  the  figures  above,  the  ideal  reflection  point  is  slightly  off  from 
the  identified  peak  of  each  edge  reflection.  The  actual  edge  reflection  times  for  the  6” 
plate  and  the  4”  flange  plate  are  .86ns  and  .57ns  respectively  for  the  thick  sample  and 
.83ns  and  .56ns  respectively  for  the  thin  sample. 

4.2.2  Signal  Subtraction 

Signal  subtraction  has  the  benefit  of  isolating  the  signal  and  does  not  broaden  the 
responses.  However,  the  effectiveness  of  the  subtraction  depends  heavily  on  the  quality 
of  the  initial  approximation  of  the  0th  order  response.  For  this  research,  a  3rd  order 
polynomial  fit  was  chosen  based  upon  heuristic  examinations  of  the  frequency-domain 
data  and  how  well  the  transformed  fit  removed  the  0th  order  reflection.  The  figures  below 
present  time-domain  data  from  the  measured  samples.  Both  a  5.55mm  and  a  4.39mm 
thick  sample  of  Plexiglas  were  used.  The  frequency  ranges  were  10.0  GHz  to  16.0  GHz 
and  10.0  GHz  to  18.5  GHz  respectively.  Again,  the  thicker  sample  has  a  reduced 
frequency  range  to  avoid  exciting  the  1st  order  mode.  The  data  covers  the  6”  and  4” 
plates  of  square  and  circular  geometry.  Zero  padded  data  uses  a  1601  point  zero  pad. 
Once  again,  as  a  frame  of  reference,  the  original  data  collected  from  the  NWA  is  included 
in  each  figure. 

As  a  final  note,  the  algorithm  used  to  process  the  data  designated  the  largest 
response  remaining  after  the  signal  processing  as  the  edge  reflection.  The  only  case  in 
which  a  large  response  may  be  ruled  out  is  if  it  occurred  before  the  time  at  which  the 
reflection  could  have  happened  in  free  space.  The  reason  for  this  is  to  evaluate  the 
method  assuming  no  a  priori  knowledge  of  the  material’s  relative  complex  permittivity. 
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|S11(t)| 


Figure  4-6:  (a)  Non  zero  padded  and  (b)  zero  padded  edge  responses  for  a  6” 
circular  plate  measuring  a  5.55mm  thick  sample 


4-9 


|S11(t)|  |S11(t)| 


Time  (ns) 


(a) 


Figure  4-7:  (a)  Non  zero  padded  and  (b)  zero  padded  edge  responses  for  a  6” 
circular  plate  measuring  a  4.39mm  thick  sample 
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Figure  4-8:  (a)  Non  zero  padded  and  (b)  zero  padded  edge  responses  for  a  4” 
circular  plate  measuring  a  5.55mm  thick  sample 
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Figure  4-9:  (a)  Non  zero  padded  and  (b)  zero  padded  edge  responses  for  a  4” 
circular  plate  measuring  a  4.39mm  thick  sample 
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Figure  4-10:  (a)  Non  zero  padded  and  (b)  zero  padded  edge  responses  for  a  6” 
square  plate  measuring  a  5.55mm  thick  sample 
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Figure  4-11:  (a)  Non  zero  padded  and  (b)  zero  padded  edge  responses  for  a  6” 
square  plate  measuring  a  4.39mm  thick  sample 
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Figure  4-12:  (a)  Non  zero  padded  and  (b)  zero  padded  edge  responses  for  a  4” 
square  plate  measuring  a  5.55mm  thick  sample 
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Figure  4-13:  (a)  Non  zero  padded  and  (b)  zero  padded  edge  responses  for  a  4” 
square  plate  measuring  a  4.39mm  thick  sample 
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For  the  circular  plates,  the  signal  subtraction  method  works  very  well  in 


indentifying  the  edge  reflection.  The  thicker  samples  provide  less  accurate  data.  The 
zero  padding  provides  no  distinct  advantage  in  identifying  the  edge  reflection.  The 
method  was  not  effective  when  square  flange  plates  were  used.  As  is  seen  in  the  above 
figures,  the  square  flange  plates  had  responses  that  were  widely  dispersed  making 
accurate  identification  of  the  edge  reflection  difficult. 

The  table  below  gives  a  summary  all  of  the  reflection  times  and  amplitudes  from 
each  of  the  figures  above.  These  values  will  be  used  in  calculating  the  relative  complex 
permittivity. 

Table  4-2:  All  the  edge  response  times,  the  difference  of  the  response  times  from  the 
ideal,  and  the  amplitudes  for  each  measurement  setup  processed  using  the  signal 
subtraction  method 


Zero  Pad 

Time  (ns) 

Diff  from  Ideal  (ns) 

Amplitude 

6" 

circular 

Thick  Sample 

N 

0.831 

0.011 

0.0471 

Y 

0.851 

0.031 

0.0456 

Thin  Sample 

N 

0.821 

0.001 

0.1026 

Y 

0.849 

0.029 

0.0913 

4" 

circular 

Thick  Sample 

N 

0.499 

-0.047 

0.0406 

Y 

0.574 

0.028 

0.0459 

Thin  Sample 

N 

0.587 

0.041 

0.1133 

Y 

0.562 

0.016 

0.0989 

6"  square 

Thick  Sample 

N 

0.831 

0.011 

0.0089 

Y 

0.851 

0.031 

0.0085 

Thin  Sample 

N 

1.291 

0.471 

0.0199 

Y 

1.267 

0.447 

0.0181 

4"  square 

Thick  Sample 

N 

0.665 

0.119 

0.0100 

Y 

0.592 

0.046 

0.0122 

Thin  Sample 

N 

0.704 

0.158 

0.0243 

Y 

0.549 

0.003 

0.0221 
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4.3  Final  Extraction 


After  the  edge  reflections  have  been  identified  via  one  of  the  two  signal 
processing  methods  addressed  previously,  the  data  can  be  used  to  calculate  the  relative 
complex  permittivity  for  the  sample  material.  The  data  will  be  presented  against  industry 
standard  values  of  s'  and  s"  for  Plexiglas  which  are  defined  here  as  s'  =  2.6  and 
s"  =  0.0150  [12].  The  sample  materials  used  for  these  measurements  were  a  5.55mm  and 
a  4.39mm  thick  sample  of  Plexiglas.  Frequency  ranges  were  10.0  GHz  to  16.0  GHz  and 
10.0  GHz  to  18.5  GHz  respectively.  As  previously  mentioned,  the  thicker  sample  has  a 
reduced  frequency  range  to  avoid  exciting  the  1st  order  mode.  Both  the  6”  and  4”  square 
and  circular  plates  were  used  to  take  measurements. 

4.3.1  Simplified  Extraction 

The  equations  for  the  simplified  extraction  were  detailed  in  (3-22)  and  (3-23).  In 
calculating  the  uncertainty  for  (3-22),  using  (4-2)  yields 


where  d2w  is  the  expected  value  for  the  distance  and  t2w  is  the  expected  value  for  time. 
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4.3.1. 1  Kaiser  Windowing 

The  results  from  applying  a  simple  extraction  to  data  from  the  NWA  that  has  been 
processed  using  Kaiser  windowing  are  displayed  in  Figure  4-14  and  Figure  4-15.  It  is 
observed  that  s'  falls  within  two  standard  deviations  (i.e.,  95%  confidence  assuming  a 
Gaussian  distribution),  s"  is  within  confidence  intervals  in  the  thicker  sample;  however, 
because  the  thicker  sample  has  a  reduced  bandwidth,  the  uncertainty  for  s'  is  greater. 
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Figure  4-14:  er  of  a  5.55mm  sample  of  Plexiglas  using  the  simplified  extraction 


and  processed  with  Kaiser  windowing. 
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Figure  4-15:  er  of  a  4.39mm  sample  of  Plexiglas  using  the  simplified  extraction  and 
processed  with  Kaiser  windowing. 

4.3. 1.2  Signal  Subtraction 

With  the  edge  reflections  identified  from  data  that  has  been  processed  using  the 
signal  subtraction  method,  a  simple  extraction  can  be  performed.  The  extracted  data 
covers  both  a  5.55mm  and  a  4.39mm  thick  sample  of  Plexiglas.  Frequency  ranges  were 
10.0GHz  to  16.0GHz  and  10.0GHz  to  18.5GFIz  respectively.  The  thicker  sample 
required  a  reduced  frequency  range  in  order  to  avoid  exciting  the  1st  order  mode  in  the 
parallel-plate  region.  The  probe  setup  used  both  the  square  and  circular  geometries. 

The  extracted  results  are  presented  below. 
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Figure  4-16:  sr  for  a  5.55mm  thick  sample  of  Plexiglas  using  circular  plates  with  no 
zero  padding 
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Figure  4-17:  sr  for  a  5.55mm  thick  sample  of  Plexiglas  using  circular  plates  with 
zero  padding 
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Figure  4-18:  er  for  a  4.39mm  thick  sample  of  Plexiglas  using  circular  plates  with  no 
zero  padding 


Figure  4-19:  sr  for  a  4.39mm  thick  sample  of  Plexiglas  using  circular  plates  with 
zero  padding 
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Figure  4-20:  £,  for  a  5.55mm  thick  sample  of  Plexiglas  using  square  plates  with  no 
zero  padding 
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Figure  4-21:  e,  for  a  5.55mm  thick  sample  of  Plexiglas  using  square  plates  with  zero 
padding 
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Figure  4-22:  e,  for  a  4.39mm  thick  sample  of  Plexiglas  using  square  plates  with  no 
zero  padding 
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Figure  4-23:  £,  for  a  4.39mm  thick  sample  of  Plexiglas  using  square  plates  with  zero 
padding 
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Examining  the  figures  above,  when  using  square  flange  plates,  the  data  for  the 


real  part  of  sr  is  almost  twice  the  industry  standard.  In  some  applications,  this  may  be  an 
acceptable  tolerance;  however,  in  most,  it  is  not.  In  contrast,  the  circular  flange  plate  data 
is  all  within  95%  confidence  intervals.  In  the  best  case,  the  uncertainty  is  only  +0.28, 
which  is  considered  an  acceptable  tolerance  for  many  industry  applications. 

With  the  assumption  of  a  low  loss  material,  the  imaginary  part  of  er ,  in  contrast  to 
the  real  part,  can  be  a  multiple  of  the  industry  standard  and  still  be  considered  within 
tolerances  (dependent  on  application)  as  long  as  it  is  within  the  same  order  of  magnitude. 
From  this  perspective,  the  circular  flange  plate  measurements  appear  to  be  good. 
However,  they  are  most  accurate  using  the  thin  Plexiglas  sample  with  no  zero  padding. 
When  examining  the  averaged  s'  and  s”  ,  the  same  observations  hold  true.  This  data  is 
contained  in  Table  4-3  where  the  averaged  s'  and  s"  are  presented  along  with  the 
industry  standard  as  a  reference. 

Table  4-3:  All  averaged  z'  and  z"  data  from  each  measurement  along  with  the 
industry  standards  for  those  values. 


Sig  Process 

plate 

Sample  Size 

Zero  Pad 

Avg  e' 

Avg  e” 

Ind  Std  e' 

Ind  Std  e" 

Kaiser 

5.55mm 

N/A 

2.846 

0.021 

2.6 

0.015 

4.39mm 

N/A 

2.699 

0.032 

2.6 

0.015 

Circular 

5.55mm 

N 

2.42 

0.080 

2.6 

0.015 

Y 

2.834 

0.049 

2.6 

0.015 

4.39mm 

N 

2.805 

0.023 

2.6 

0.015 

Signal 

Y 

2.768 

0.028 

2.6 

0.015 

Subtraction 

5.55mm 

N 

3.261 

0.021 

2.6 

0.015 

Square 

Y 

2.927 

0.040 

2.6 

0.015 

4.39mm 

N 

5.381 

0.001 

2.6 

0.015 

Y 

4.414 

0.000 

2.6 

0.015 
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4.4  Summary 

This  Chapter  has  presented  the  results  of  using  a  simple  time  domain  based 
extraction  of  complex  permittivity  that  utilizes  the  edge  reflections  from  waveguide 
probes.  Starting  with  an  examination  of  the  uncertainty  associated  with  the  system,  the 
signal  processing  methods  were  discussed,  and  the  results  of  the  processing  were  used  in 
extracting  the  permittivity.  Two  samples  were  used  for  the  measurements,  a  5.55mm 
thick  and  a  4.39mm  thick  sample  of  Plexiglas.  Furthermore,  the  frequency  range  covered 
for  the  measurements  was  10.0GHz  to  16.0GHz  and  10.0GHz  to  18.5GHz  respectively. 
Two  types  of  probes  were  used  for  taking  the  measurements.  The  first  type  was  a  Ku- 
band  waveguide  with  either  a  6”  or  4”  square  flange  plate  attached.  The  second  type  was 
also  a  Ku-band  waveguide  with  either  a  6”  or  4”  circular  flange  plate  attached.  It  was 
found  zero-padding  did  not  affect  the  final  results  when  using  the  signal  subtraction 
method.  Furthermore,  the  results  when  using  circular  flange  plates  are  considered  to  be 
acceptable  for  use  in  most  industry  applications.  Finally,  higher  order  modes  should  be 
avoided  so  that  the  8.5GHz  bandwidth  of  the  Ku-band  waveguide  can  be  used.  Since 
temporal  resolution  was  the  primary  driver  of  uncertainty,  avoiding  a  reduced  bandwidth 
in  the  Ku-band  waveguide  will  lead  to  more  accurate  measurements  with  the  desired 
confidence  intervals. 
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5.  Conclusions  and  Recommendations 


A  simple  method  has  been  developed  that  can  produce  reasonable  results  for 
extracting  relative  complex  permittivity  from  a  low-loss,  non-magnetic,  non-dispersive 
dielectric.  Plexiglas  was  used  in  this  research  as  a  proof-of-concept,  but  the  method  can 
be  applied  to  other  materials  within  this  classification.  The  main  advantage  of  the 
method  is  that  it  is  both  non-destructive  and  computationally  simple.  This  is  of  huge 
benefit  over  other  non-destructive  methods  that  often  use  a  Method  of  Moments  solution 
which  can  be  time  consuming  and  computationally  demanding.  Furthermore,  since  the 
method  is  non-destructive,  no  time  is  lost  in  sample  preparation,  and  no  errors  are 
introduced  by  a  poorly  fitted  sample. 

5.1  Conclusions 

During  the  development  of  this  method,  two  measurement  setups  and  two  signal 
processing  methods  were  explored  as  a  means  of  accomplishing  the  measurement  and 
analyzing  the  signal.  This  was  done  intentionally  so  that  the  method  could  be  further 
refined  in  later  work  without  needing  determinations  to  be  made  regarding  plate 
geometries  or  signal  processing. 

5.1.1  Preferred  Plate  Geometry 

Based  on  the  results  presented  Chapter  4,  plate  geometry  appears  to  have  a 
profound  effect  on  the  real  part  of  the  relative  complex  permittivity.  This  effect  is  caused 
because  the  1st  order  reflection  may  no  longer  be  considered  the  edge  reflection  in  square 
flange  plates.  While  2nd  order  reflections  at  about  the  appropriate  times  are  noticed, 
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positive  identification  cannot  be  made  since  the  permittivity  is  assumed  to  be  unknown. 
Therefore  the  circular  plates  are  the  preferred  geometry  for  use  with  this  method. 

5.1.2  Preferred  Signal  Analysis  Method 

Unlike  the  plate  geometry,  the  final  extracted  data  does  not  unequivocally  show 
one  signal  processing  method  to  be  superior  to  another.  The  average  for  s'  for  Kaiser 
windowing  (with  zero  padding)  is  higher  than  the  averages  obtained  via  signal 
subtraction,  but  when  comparing  Figure  4-14  to  Figure  4-16,  the  values  for  s'  are  within 
95%  confidence  intervals.  However,  waveform  spreading  occurs  when  using  Kaiser 
windowing,  which  subsequently  diminishes  temporal  resolution.  Therefore,  Signal 
Subtraction  is  identified  as  the  preferred  signal  processing  method.  Furthermore,  the 
NWA  used  in  this  research  was  able  to  apply  the  Kaiser  window  which  added  to  the 
simplicity  of  the  method.  Other  NWA’s  may  not  possess  this  capability  which  removes 
the  key  benefit  of  using  that  signal  processing  method. 

5.1.3  Confidence  Interval 

The  temporal  resolution  uncertainty  was  based  on  the  assumption  that  there  was  a 
uniform  probability  for  the  1st  order  reflection  to  occur  somewhere  within  an  interval  that 
is  the  size  of  the  resolution  and  centered  on  the  point  indentified  as  the  measured  edge 
response.  This  is  a  valid  assumption  if  is  unknown  for  the  sample  material.  Therefore, 
the  exact  location  of  the  edge  reflection  remains  unknown  and,  subsequently,  has  an 
equal  chance  of  being  anywhere  within  the  region  identified  with  the  edge  response.  This 
leads  to  an  uncertainty  in  the  s'  measurement  of  ±0.28  when  using  a  95%  confidence 
interval  and  the  full  frequency  range  of  the  Ku-band  waveguide.  However,  the  PDF  can 
be  altered  from  a  unit  step  function  to  another  form  that  may  better  describe  the  material 
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if  there  is  some  a  priori  knowledge  of  the  material.  With  an  altered  PDF,  the  uncertainty 
can  be  reduced  leading  to  increased  precision. 

5.1.4  Final  Recommendation 

Based  upon  the  recommendations  made  above,  the  method  described  in  this 
research  can  be  optimized  in  the  following  ways.  First,  circular  flange  plates  should  be 
used  to  reduce  ambiguity  associated  with  identifying  the  edge  reflection.  Additionally, 
since  temporal  resolution  significantly  affects  the  uncertainty  in  the  measurement, 
anything  that  will  reduce  the  bandwidth  of  the  system  (such  as  having  to  lower  the  upper 
limit  of  the  frequency  range  in  order  to  avoid  generating  high  order  modes)  should  be 
avoided.  Finally,  the  signal  subtraction  method  of  signal  processing  should  be  used  to 
avoid  waveform  broadening  and  loss  of  temporal  resolution. 

However,  even  with  this  setup,  the  method  is  not  without  drawbacks.  Because  of 
the  inherent  simplicity  of  the  system,  material  type  is  limited  to  a  certain  form  of 
dielectrics.  Furthermore,  uncertainty  is  introduced  into  the  measurement  via  the  signal 
processing  and  the  assumptions  made  in  order  to  keep  the  calculations  computationally 
simple. 

5.2  Recommendations  for  Future  Research 

Research  should  be  conducted  into  increasing  the  bandwidth  of  the  measurement 
apparatus.  With  increased  bandwidth,  the  time  resolution  would  be  improved  which 
would  lead  to  more  precise  and  accurate  measurements.  Another  area  for  future  research 
is  to  looking  into  somewhat  more  rigorous  mathematical  models.  As  another  research 
area,  the  signal  subtraction  method,  which  relies  on  a  simple  polynomial  fit,  should 
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explore  the  use  of  Hermite  or  Legendre  polynomials.  Additional  signal  processing 
research  should  also  focus  on  the  use  of  other  windows  as  effective  signal  processing 
methods.  The  Kaiser  window  was  considered  to  be  representative  of  similar  windows 
(Hann,  Hamming,  etc),  but  without  a  formal  study,  the  other  windows  cannot  be  ruled  out 
as  effective  signal  processing  methods.  A  final  area  to  investigate  is  expanding  or 
modifying  the  method  to  include  magnetic  materials. 
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Appendix  A:  Acronyms 


Bandwidth  BW 

Discrete  Fourier  Transform  DFT 

Electromagnetic  EM 

Fourier  Transform  FT 

Gigahertz  GHz 

Inverse  Fast  Fourier  Transform  IFFT 

Network  Analyzer  NWA 

Nicholson  Ross  Weir  NRW 

Nano-second  ns 

Perfect  Electrical  Conducting/Perfect  PEC 

Electric  Conductor 

Probability  Density  Function  PDF 

Short-Open-Foad-Thru  SOFT 

Thru-Fine-Reflect  TRF 
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Appendix  B:  Plate  Drawings 


This  section  includes  the  drawings  used  by  the  AFIT  machine  shop  to  fabricate  the  plates  that 
were  used  to  form  the  waveguide  probes,  they  are  included  here  as  a  reference  in  case  plates 
need  to  be  fabricated  again  in  the  future. 
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Appendix  C:  Waveguide  Bandwidth 

A  Ku-band  waveguide  was  used  as  part  of  the  probe  setup.  The  operating  bandwidth 
used  was  10.0GHz  to  18.5GHz.  That  bandwidth  was  determined  based  upon  the  following 
calculations. 

Whereas  the  parallel  plate  system  relied  on  operating  in  the  0th  order  mode,  a  wave  guide 
requires  at  least  the  1st  order  mode  to  be  able  to  excite  waves.  For  purposes  of  this  application, 
the  TEjo  mode  was  used.  Therefore,  the  bandwidth  of  the  system  is  dependent  on  two  factors. 
The  first  is  TEio  mode  cutoff.  Frequencies  below  this  threshold  will  not  excite  waves  inside  the 
waveguide  and  will,  therefore,  be  useless.  This  defines  the  low  point  for  the  bandwidth.  The 
high  point  is  defined  by  the  TE20  cutoff  frequency.  This  point  is  chosen  so  that  no  higher  order 
modes  are  excited.  The  formula  for  the  cutoff  frequency  for  a  TEmn  mode  of  a  freespace  filled 
waveguide  is  given  by  Balanis  [4]  as 


fc 

2  n 


f  \  2 

'  mn  ' 


/  \  2 

'  H7T  ' 


V  a  J 
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where  a  and  b  represent  the  width  and  height  of  the  aperture  waveguide  respectively.  For  a  Ku- 
band  waveguide,  those  dimensions  are  0.31 1”  and  0.622”. 

When  calculated,  the  cutoff  frequency  for  the  TE10  mode  is  9.49GHz,  and  the  cutoff 
frequency  for  the  TE20  mode  is  18.975GHz.  Typically,  the  lower  bound  is  multiplied  by  1.25 
and  the  upper  bound  by  0.95  in  order  to  ensure  that  the  waveguide  is  operating  only  in  the  range 
of  the  desired  mode.  For  a  Ku-Band  waveguide,  this  is  usually  defined  as  approximately  12GHz 
to  18GHz.  However,  since  bandwidth  is  vitally  important  to  this  method,  a  compromise  was 
made  in  order  to  expand  the  bandwidth  beyond  the  6GHz  covered  by  the  traditionally  defined 
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range.  This  compromise  settled  on  using  a  range  from  10GHz  to  18.5GHz.  Theoretically,  the 
boundaries  could  be  expanded  to  cover  9.5GHz  to  18.9GHz;  however,  this  would  introduce  a 
risk  of  higher  order  modes  corrupting  the  data. 
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